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The eddy viscosity for a turbulent compressible fluid with a relativistic equation of state is derived. 
Compressibility allows for sound modes, but the eddy viscosity in the shear mode is found to be 
the same as for incompressible fluids. For two space dimensions (which is the relevant case for 
the dynamics of relativistic heavy-ion collisions) the eddy viscosity in the shear mode is negative, 
reducing the effective viscosity below its microscopic value. This could explain the tiny viscosity 
found at RIflC. Implications for the experimentally accessible elliptic flow coefficient at the LHC 
are speculated on. 



Efforts to understand the bulk physics of the ongoing 
experimental program at the Relativistic Heavy-Ion Col- 
lider (RHIC) have led to an interest in the theory of rela- 
tivistic fluid dynamics, both with and without shear vis- 
cosity. The curious fact that ideal fluid dynamics is able 
to describe the experimentally measured particle spectra 
amazingly well [J, 0, 0, has led to the hypothesis 
that the shear viscosity 77, or rather the relevant dimen- 
sionless quantity rj/s involving the entropy density s, has 
to be extremely small [1, Indeed, values of 77/s ^ 1 
(in natural units kB = fi = c = I) obtained by solving 
Quantum-Chromodynamics (QCD) in the weak coupling 
expansion were seemingly ruled out 9]. Interest in 
RHIC spread much beyond nuclear physics when it was 
realized that a calculation of 77/s by string theory meth- 
ods for a thermal gauge theory in the strong coupling 
limit [lOj gave 77/s — about an order of magnitude 
smaller than the weak coupling QCD result. 

While this result had originally been derived not for 
QCD but for the J\f — 4 supersymmetric Yang-Mills the- 
ory, it was later conjectured that all relativistic quan- 
tum field theories at finite temperature and zero chemical 
potential obey 77/s > No laboratory fluid seemed to 
violate this bound until recently, when it was shown [l^ 
that the experimental data from top RHIC energies ac- 
tually favored 7 7/s ~ (another group argues for even 
smaller values [l3|). While addressing the caveats out- 
lined in 12] will ultimately allow to decide whether RHIC 
violates the bound 7;/s > ^ or not, it is probably fair to 
conclude that 77/ s at RHIC is apparently extremely small 
(see also [i3, El, El)- 

The Reynolds number for a relativistic fluid at tem- 
perature T can be estimated as Re ^ ffTL, where L is a 
typical length scale. For gold collisions at RHIC, taking 
L to be the radius of a gold nucleus L 6 fm and the 
QCD scale as temperature (T ~ 200 MeV) imply 
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where the value for 77/s from |12l | was used. One is thus 
led to wonder about the possibility and consequence of 
fluid turbulence in heavy-ion collisions at RHIC. 



The natural starting point are the relativistic fluid dy- 
namic equations in the presence of shear viscosity 

{e+p)Du^' = V>-A^D^H"'3, 

De = -(e+p)V^7^^ + iH^''(V,u^}, (2) 

where e,p are the energy density and pressure, is the 
fluid four velocity obeying m'^m^ = 1 and H^'' is the shear 
tensor which to lowest order in gradients [s^] is H'^'^ = 
ri{W^u'^). The remaining symbols in Eq. ([2]) are: Da, the 
covariant derivative, D = u^^D^j,, = AJ^D^, A'"'' = 
gf"" - u>'u'' and {Vu") = Vt^u" + Vu'' - fA'^'^Vaw". 
The metric signature is g'^" — (+, — , — ) and (+, — , — , — ) 
for d = 2 and d — 3 space dimensions, respectively. The 
equation of state p — p[e) is used to close the system ([2]). 
In what follows, an equation of state with a constant 
speed of sound squared = dp/de will be considered. 

Eq. ^ is much more non-linear than the non- 
relativistic Navier-Stokes equation for incompressible flu- 
ids, which is the canonical starting point for turbulence 
theory. In order to make progress, it is therefore useful 
to consider the approximation where both 77/s and the 
fluid three velocity are small, 77/s zi* ~ C>{o) ^ 1. 
Introducing the Fourier transform for d + 1 dimensional 
space-time 

lne(t,x)= J j^--^^exp[-~iujt + tk-x\lne{LJ,k) (3) 
and similarly for u*, one flnds (in flat space) from Eq. ([2]): 
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where K = (u;,q), and ^{/{e + p) was approximated by 
its space-time average, v = Jdtd'x.r]/{e + p). This ap- 
proximation ~ though not strictly justified - will allow 
to connect to literature on non-relativistic turbulence, 
where v is referred to as kinematic viscosity. 

Following the Yakhot-Orszag approach to turbulence 
[3], Eq. (HI is replaced by a more general equation 

{D^\K)f < = f{K) f u^u'k^qP^'^^{K, Q), 

^ Q 

(5) 

where p is a random force used to model turbulent stir- 
ring at high wavenumbers fc > A. Its correlator in d space 
dimensions is taken to be Gaussian 

< P{K)P{K') 2ak-'^{2T:Y+^5'^5{K + K') (6) 

with strength a and zero mean < P{K) >= 0. Sphtting 
= u'^ + into a low and high wavenumber part 
(u'^< = u^6i(A-fc)), the aim is then to derive an averaged 
equation for the geometric flow u'^^ in the presence of the 
fiuctuating 



-W{K)0{K-k) / P"™(i^,g)x 



( 



m,< /.< 



'■K-C 



. m.> /,> ^ 



(7) 



To calculate the mean < u™' u^_q >, one solves Eq. ([5]) 
for recursively in [18] . Using Eq. ^ and focusing 
on the term linear in it follows that 
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Defining the projectors A% = S'^ - B% and B%^!^, 
one can decompose the propagator Dq into a shear and 
a sound mode (similar to the Kovasznay modes [l^). 
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where J^ ^ = /g 6'(g - A)6'(|k - q| - A). The fre- 
quency integrations are readily performed and in the 
small wavenumber limit a;, k — > one finds 
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where only terms larger than 0(1^) ~ C(o) were kept. 
Again, i?" can be decomposed as i?" = r^A^ -t-rsB^f . 
From Eq. pUj) this implies that 
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where D*-' — {Dq ) ^ (if) + R^^ is the propagator involv- 
ing turbulence corrections, 

(Z?«^ )-i = (a^ + + (I3j, + rB)B%. (13) 

To evaluate in the small k limit, one contracts R^°-A''^ 
and upon shifting q ^ q -I- ^k and using the symmetries 
of the integral obtains (see also [3]) 
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Interestingly, this result is identical to that for incom- 
pressible fluids [Til and modifications for compressible 
fiuids only arise in the sound channel of the propagator 
(which is absent for incompressible fluids). The equation 
for the full propagator Eq. (fT3|) then implies that 
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where the effective turbulent viscosity Uf,f / — v ~\- Vf,, 
being the sum of microscopic and eddy viscosity 



The inversion of Dq is then straightforward and one 
finds 
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has been introduced, v^ddy refiects the effect of turbu- 
lence acting on length scales smaller than A~^ onto the 
geometric flow at scales larger than A~^. Curiously, note 
that Veddy IS positive for d = 3 space dimensions, while 
it is negative for d = 2. In other words, the effective 
viscosity in a system exhibiting two-dimensional turbu- 
lence is smaller than the microscopic viscosity, due to 
the presence of a negative eddy viscosity. This interest- 
ing phenomenon has been known for decades, starting 
with the work of Starr and Kraichnan 2^ 2l[ (see also 
2^ . 2^ for other approaches). 



In heavy-ion collisions, one typically follows Bjorken 
in assuming boost-invariance in the longitudinal direc- 
tion [i^l- Therefore, the dynamics of relevance for fluid 
dynamics is effectively two-dimensional. Hence, if fluid 
turbulence develops, one expects the eddy viscosity to be 
negative and thus the fluid at RHIC would appear more 
ideal than it is based on its microscopic viscosity. 

The physics of the appearance of a negative eddy vis- 
cosity seem to be tied to the phenomenon of inverse en- 
ergy cascade. While in three dimensional turbulence, 
energy generally cascades down to smaller and smaller 
length scales until it is finally dissipated into heat (reg- 
ular cascade), in two dimension this process can seem- 
ingly reverse: energy is transferred from smaller to larger 
length scales (see also the discussion in Ref. [llj). A key 
ingredient seems to be that vorticity is conserved in two- 
dimensional incompressible fluids, allowing small eddies 
to convey energy to larger eddies. For compressible rel- 
ativistic fluids, vorticity is no longer conserved and 
hence one can expect the inverse cascade to "leak" en- 
ergy. Indeed, working out the turbulent correction to the 
sound mode of the propagator, one finds an which is 
identical to except that in (fT4|l one has to replace 
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This implies that > for both d = 2,3 and ^ ^. 

It seems that in two-dimensional compressible fluid 
turbulence the shear mode gets less dissipative while for 
sound the converse is true. This means that it could 
be difhcult to decide whether to expect an increase or 
decrease of effective viscosity in different experimental 
observables. The very small apparent viscosity at RHIC 
has been extracted mainly from the experimental observ- 
able called "elliptic flow" 12, 2^. Invoking Morkovin's 
Hypothesis that under certain conditions in compressible 
turbulence "the essential dynamics of these shear flows 
will follow the incompressible pattern" il9-], in the fol- 
lowing it will be speculated that elliptic flow is mostly 
influenced by the negative eddy viscosity in the shear 
mode. Then, is it possible to verify or falsify the idea of 
fluid turbulence in relativistic heavy-ion collisions? 

Clearly, some other experimentally accessible predic- 
tion from Eq. p6)) is needed. To do so, first introduce 
the dimensionless coupling parameter A = -jjsjyc so that 
^e// = Mtt)' "^hsre d = 2 has been used in Eq. p6)l . 
Now a renormalization group improvement is performed 
by solving the set of equations [18| 
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where to leading order in A, v has been replaced by Veff 
and a was assumed to be independent of A [l^ . Using 

1/3 

VeSf\\=Q = V, the result is Veff = u {l - ^) 
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FIG. 1: Effective versus standard viscosity for A = 10~^. De- 
creasing 77/s from large values, the effective viscosity changes 
rapidly to extremely small values close to the critical viscosity 
rj/s = (A)^''^. Decreasing 77/s further may even result in nega- 
tive effective viscosity. However, at this point the calculation 
breaks down and the result cannot be trusted. 



other words, one expects the effective ratio 77/s to behave 
as 
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where for convenience A = 
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has been introduced. 
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It may be that A still depends on -q/ s: for the following 
qualitative discussion, this should not matter unless A is 
proportional to 77/s with a power greater or equal three. 
A will typically be rather small: since a has to be propor- 
tional to the temperature (the only other dimensionful 
scale except A), one expects A ^ (T/A)**. Assuming A 
corresponds to the smallest resolved length scale of cur- 
rent heavy-ion collision simulations [l3 | one has A 1 
GeV, so using again T - 200 MeV gives A - lO"'''. 

The qualitative behavior of the effective viscosity is 
shown in Fig.[l]for A = 10~^. It can be seen that close to 
the critical viscosity the effective viscosity becomes very 
small and there is some indication that even slightly neg- 
ative values of effective viscosities can be achieved. It 
is expected that higher non-linearities subsequently ren- 
der the effective viscosity again positive [23[ . Since these 
terms have been ignored in the above calculation, the 
result in Eq. cannot be trusted below the critical 
viscosity. Nevertheless, the strong suppression of the ef- 
fective viscosity close to the critical viscosity suggested 
by Eq. (fT8|) should be experimentally accessible, e.g. by 
a measurement of the elliptic flow at the upcoming Large 
Hadron Collider (LHC). The magnitude of elliptic flow 
is known to decrease for increasing 77/s Turning the 
argument around, if fluid turbulence is in operation at 
RHIC, and elliptic flow is indeed mostly affected by the 
negative eddy viscosity in the shear mode, the effect at 
the LHC should be even more pronounced: from Fig[T] 
one would expect the elliptic flow to increase beyond the 



4 



RHIC values, maybe even beyond the "ideal hydrody- 
namic limit", if a negative effective viscosity is realized 
for an extended period of time. Note that more conven- 
tional approaches [2^, 27 1 call for a decrease of elliptic 



flow at the LHC, while extrapolations from existing data 
do indicate an increase 
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It should be stressed that the above calculation of the 
eddy viscosity for fluid turbulence in relativistic heavy- 
ion collisions can be regarded as qualitative at best. How- 
ever, given interest in the community, there are many 
ways to improve or strengthen Eq. p^ . e.g. by relax- 
ing various assumptions made in the derivation or build- 
ing upon the considerable knowledge from studying two- 
dimensional incompressible fluid turbulence (see e.g. [2^ 
and references therein). 

Finally, it should be pointed out that the concept of 
anomalous turbulent viscosity in the context of heavy- 
ion coUisions has been already suggested in Ref. [30|. 
However, the discussion in [30| is based on assuming 
the presence of plasma turbulence, which is somewhat 
different than the fluid turbulence picture outline here. 
(Non-Abelian) Plasma turbulence seems to occur as a 
consequence of plasma instabilities 31, H, 33, 34|. How- 
ever, given current estimates of initial plasma parame- 
ters it seems somewhat unlikely that plasma instabilities 
set in early enough 34, 3^ for plasma turbulence to be 
relevant at RHIC. Nevertheless, the spectrum of initial 
fluctuations studied in this context [36| could also be of 
relevance for fluid turbulence. 

To summarize, following the Yakhot-Orszag approach 
to turbulence the eddy viscosity in a compressible fluid 
with a relativistic equation of state was calculated. The 
result differs from incompressible fluids by the presence 
of a sound channel in the fluid dynamic propagator. In 
the sound channel, the eddy viscosity is positive for both 
two and three dimensions. In the shear channel, the eddy 
viscosity is found to be identical to that of incompressible 
fluids, and happens to be negative for two space dimen- 
sions. Given the Reynolds number for RHIC experiments 
is expected to be much larger than unity, and that the 
relevant fluid dynamics is essentially two-dimensional, it 
may be that this explains the apparently tiny viscosity 
over entropy density ratio extracted from experimental 
data. If this is the case, the model would predict a sub- 
stantially increased elliptic flow at the LHC, which is 
experimentally falsiflable in the near future. 
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